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\Q ABSTRACT 

o 

Aims. We analyze the impact on the Galactic nitrogen abundances of using a new set of low and intermediate mass star yields. These yields 
' have a significant yield of primary nitrogen from intermediate mass stars. 
II Methods. We use these yields as an input to a Galactic Chemical Evolution model and study the nitrogen abundances in the halo and in the 
© ' disc, and compare them with models obtained using other yield sets and with a large amount of observational data. 

Results. We find that, using these new yields, our model adequately reproduce the observed trends. In particular, these yields solve the 
C/2 . historical problem of the evolution of nitrogen, giving the right level of relative abundance N/O by the production of a primary component 
in intermediate mass stars. Moreover, using different evolutionary rates in each radial region of the Galaxy, we may explain the observed N 
dispersion. 
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1. Introduction 

Most elements are creat ed in the interiors of s tars by nucle- 
osynthesis processes (see lWallerstein et al J 19971 for a review), 
starting with hydrogen and progressing toward heavy elements. 
These processes are called primary production. Some elements, 
however, can be formed from nuclei heavier than hydrogen 
originally present in the star. They are called secondary. This 
is the case of nitrogen, that can be created during the CNO 
cycle using seeds of original carbon and/or oxygen. From a 
theoretical point of view, it has b een considered t hat massive 
stars produce secondary nitrogen (Peimbertl ll987h . while low 
and intermediate mass (LIM) stars have mechanisms, like the 
third dredge-up and the Hot Bottom Burnin g processes, to pro- 
duce both, primary and secondary nitrogen ( Ed munds & P agel 
ll978UAlloin et alJl979h . The third dredg e-up event is a conse- 
quence of the thermal pulses in the star, and transport C and He 
to the outer layers. The Hot Bottom Burning occurs when the 
CNO cycle takes place at the base of the convective envelope. 

Send offprint requests to: Marta Gavilan 



Observationally, there are several open questions about the 
primary or secondary character of nitrogen that up to now re- 
main unsolved. When N and O data are represented as log(N/0) 
vs log(0/H), including the galactic stars, Hn re gions for the 
Milky Way Galaxy (MWG), external galaxies (iGarnett et a"fl 
1 19951 Il999t Ivan Zee et alJll998at llzotov & Thuanlll999l) and 
the high redshift data (|Pettini et alJ 120021 [Prochaska et alJ 
l2002t ICenturion et alJl2003l and references therein), a clear 
positive slope appears for abundances larger than 12 + 
log(0/H) = 7.8 - 8 dex which indicates a secondary behavior, 
but the plot shows a flat slope for low metallicities that can only 
be explained with a primary component of nitrogen. Taking 
into account that this flat slope occurs for low abundance s, the 
first i dea prop osed , shared by some auth ors (Page l et. all 19791 
biaz & Tosilll98d bahmen et alJll995tl . is that observations 
would be reproduced if the nitrogen ejected by massive stars 
would be primary, while intermediate mass stars might have 
both primary and secondary components. 

Thus, some authors have tried to look for mechanisms 
that explain how massive stars could produce primary nitro- 
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gen. This is the case of lMevnet&MaedeJ J2002h that have 
recently proposed rotation as a possible source of primary 
nitrogen, since low metalli city stars show a bigg er rotation 
than high metallicity ones. Chiappini et al. (2003a) have used 
these yields in their chemical evolution models, concluding 
that they are only a lower limit for the primary nitrogen pro- 
duction since the Hot Bott om Burning i s not considered in 
their calculation. In fact, Chiap pini et alJ ll2005h find that an 
extra- production of N in low metallicity massive stars by 
a large factor, between 40 and 200 along the mass range, is 
necessary to explain the data of very metal-poor halo stars 
since these yields do not produce a sufficient amount of pri- 
mary N. Moreover, if the production of primary nitrogen would 
proceed from massive stars, the left side of the (N/O) vs 
(O/H) plot should n ot show any scatter. Although some authors 
claim to observe dlzotov & Thuanlll99 9: Pilvu gin et alJ l2003) 
this no-scatter, recent observations from low metallicit y ob- 
jects JPettini et al.l2002llProchaska et al.l2002HCenfuri6n etail 
l2003t llsraelian et al.ll2004 LSpite et al.ll2005h do show a clear 
dispersion. 
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ISerrano & Peirribertl lll983l) already claimed that a sec- 
ondary production by intermediate mass stars must exist and 
suggested that the zero slope may be explained by two factors: 
1) a delay in the ejection of N to the ISM due to the different 
mean-lifetimes of stars and 2) the gas infall effects. The advan- 
tage of taking a delay into account is that the great data scatter 
can be explained by considering different evolutionary states 
for each galaxy and, therefore, this possib ility has been sup- 
ported by a large number o f authors dWa Costas & Edmunds 
1993 | IPilvuginll 19921 Il993l Ivan Zee et alJll998bt iHenrv et alJ 
2000). These last ones also include gas flows -infall or/and 
outflow-, and low efficiency for the star formation rate (the 
equivalent mechanism to produce a delay) in the low evolved 
regions, in order to reproduce the flat slope in the (N/O) vs 
(O/H) plot. They conclude that the secondary production of ni- 
trogen should dominate in high metallicity environments while 
the primary one should act at low metallicities. 

Some new yiel ds for LIM stars have been given in 
hereinafter Paper I), where they were ad- 
equately evaluated and calibrated by using them in a Galaxy 
chemical evolution model. It was shown that the results about C 
and O abundances adequately reproduce the Galactic and Solar 
Neighborhood data. The purpose of this work is to analyze the 
impact of these stellar yields on the nitrogen abundances. In 
particular, we check, using the same prescriptions of Paper I, 
if the contribution to N given by these yields for LIM stars is 
sufficient to justify the amount of primary nitrogen the obser- 
vations point out. 

We describe the yields in Section 2, analyzing in particular 
the primary and secondary components of the nitrogen produc- 
tion. In section 3 we describe briefly the chemical evolution 
model. Section 4 is devoted to the results, and the conclusions 
are presented in section 5. 



Fig. 1. Total yields (BU) of: a) n C and b) U N produced by 
LIM stars for different metallicities following label on panel 
b), expressed as log(Z/Z ). 

2. Low and intermediate mass yields: The 

secondary and primary components of nitrogen 

The aim of this work is the study of the nitrogen behavior, using 
the same set of yields as in Paper I, that we call BU yields. For 
comparison purposes we also take th e LI M stars yields from 
Ivan den Hoek & Groenewegenl dl997l) and lMarigold200ll) that 
we call VG and MA, respectively. The complete table of yields 
BU was already given in Paper I for five metallicities: -0.2, - 
0.1, 0.0, +0.1 and +0.2, e xpressed as log(Z/Zp), wh ere solar 
abundances are taken from lGrevesse & Sauvall (11998) l . 

We summarize the behavior of the carbon and nitrogen 
yields for LIM stars, as shown in Fig. [0 In panel a) we see 
that 12 C yield is extremely small for stars with mass lower 
than 2 M Q , since they do not experience third dredge up events. 
However, stars begin to suffer these kinds of events for smaller 
masses at lower metallicity. In other words, in the low mass 
range, the metallicity and the 12 C yield are anti-correlated. 
When the stars have enough mass to undergo Hot Bottom 
Burning (HBB), the 12 C yield drops abruptly because of the 
conversion of carbon into nitrogen. The l4 N yield presents a 
local maximum in the mass range from 3.5 to 5 M Q , depending 
on the metallicity, then decreases before beginning to increase 
again as a function of stellar mass. The largest amount of ni- 
trogen is produced by stars of intermediate mass because HBB 
and the 2 nd dredge-up occur only in stars with M > 3.5 - 5M . 

1 The use of these solar abundances implies that Z Q = 0.02. 
Recently, lAsplund et alJ 120051) have obtained lower abundances 
which lead to a value Z Q = 0.012. However, these new determi- 
nations are still questioned by some au thors iBa hcall et al J 12005: 
iDrake & Testall2005l:lAntia & BasJEo 05 ) because they do not fit the 
helioseismological constraints. 
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Fig. 2. The ratio of the yield of primary N to the total yield 
of 14 N as derived from different author sets: BU, VG and MA. 
Panel a) a comparison of the solar metallicity yields for the 
three sets as labeled. Models for two different compositions, as 
labeled, are represented in panels b), c) and d) for BU, MA and 
VG, respectively. 

As the HBB increases the luminosity and the mass-loss rate, 
stars that suffer this process have shorter TP-AGB lifetimes. 
The local maximum occurs in the transition between stars with 
HBB and those without. The increase at higher masses is due to 
the shortened time between third dredge-up events. The yields 
at the lowest masses are due to the 1 " dredge-up. 

The most important difference among the used yields re- 
sides in the contribution of primary and secondary components 
of nitrogen by LIM stars. In Fig.[5]we represent the fraction of 
primary l4 N for the three used sets, as labeled, as a function 
of mass (M < 8 M , except for MA for which M < 5 M ). In 
panel a) we show the results for solar abundances. All of them 
show a similar behavior with a maximum for masses around 
3.5-4 M . 

We must clear some points about the components of N. The 
only difference between primary and secondary nitrogen is the 
origin of the carbon atom producing it. Although the idea is 
conceptually clear, it is not so simple to separately compute 
both components. Thus, although BU and MA give the two 
components separately for each model, VG do not. These au- 
thors, however, show their yields in each phase of stellar evo- 
lution. If we consider that all the nitrogen created in the AGB 
phase is primary, about ~ 90 % of the N ejected by LIM stars 
will be primary. This is sometimes assumed when these yields 
are used. This hypothesis, that we call AGB technique leads to 
a primary N component excessively large and is not totally ad- 
equate. 

Let us return to the definition: secondary N proceeds from 
the burning of original 12 C. If a fraction of the original car- 
bon is burned in the pre-AGB phase, it produces secondary 
N. Sometimes, this gives a negative 12 C yield. But, not all 
the initial carbon is consumed before the AGB phase. If we 
take, as an example, a star of 4M of solar abundance, that is 
with X( l2 C) = 0.28 x 10~ 2 , it has an initial 12 C abundance 



4M X( 12 C) = 1. 12 xlO" 2 M The pre-AGB phase carbon yield 
is yC\2 pre = 0.300 x 10~ 4 , so the mass of this element present 
in the star before the AGB begins is: 

M( 12 C) = yC\2 pre M ini + M md X( n C) (1) 

where M en d is the mass of the star at the end of this first 
phase: 3.95 M . Therefore, there is a mass ~ 1.118 x 1O~ 2 M 
of 12 C, from which M end X{ n C) ~ 1. 1 02 x 1O~ 2 M corresponds 
to original carbon. This implies that a quantity of the origi- 
nal carbon is still available to form nitrogen in the following 
phases. Thus, a fraction of the total nitrogen produced in the 
AGB phase (given by the addition of the two values given by 
VG in their tables denoted AGB yields and final AGB yields) 
may be secondary. In order to calculate this component from 
the total AGB yields we use the fraction, called r, between the 
secondary to the total nitrogen yield, r = 14 Ns/ l4 N. Taking 
into account that the secondary N proceeds from the existing 
carbon used as a seed, we assume that r is equal to the ratio 
between the old carbon and the new plus old carbon: 

= U N S _ l2 C oU 
r l4 N 12 C„ W + 12 C new 

This method (hereinafter called r method), may only be ap- 
plied to stars which suffer the HBB and produce primary lA N, 
that is, those for which the core mass before the HBB is larger 
than M HBB - O.8M , usually stars with M > 3.5 - 4M . 
Otherwise, the nitrogen yield is all secondary. The results of 
VG shown in Fig.|3proceed from this calculation. 

We have then computed the integrated yields of l4 N pro- 
duced by LIM stars that we present in Fig. [3] as a function of 
metallicity Z. In panel a), we represent the BU results as solid 
circles to which we have performed a least-squares fit shown 
by the solid (red) line. This integrated yield for l4 N, equivalent 
to the yield produced by a single stellar population, is located 
between the two other sets in this panel, with a similar depen- 
dence on Z that VG 2 but with lower absolute values. 

More significant however, is the dependence on metallic- 
ity of the ratio of primary to total integrated yields, l4 Np/ l4 N, 
shown in panel b). This ratio increases for decreasing metal- 
licity for all sets, as expected, although for VG the integrated 
yield is quite different if we consider AGB technique than r 
method. 

This metallicity effect can be easily explained: low metal- 
licity stars have smaller radii and take longer to reach super- 
winds, so they have more time to experience more third dredge- 
up events than solar metallicity stars. As a consequence, they 
have more fresh 12 C in their envelopes and they can make more 
primary nitrogen by the HBB process. On the other hand, due 
to the lower amount of original carbon, they produce less sec- 
ondary nitrogen. For VG yields the ratio is almost constant at 
a value of 15% when the r technique, represented by the short- 
dashed (blue) line, is used, for metallicities greater then 0.01, 
although it also increases for metallicities lower than this value. 

2 The integrated yields for VG are slightly different than those ob- 
tained bv lHenrv et alj l200Ct) due to the Initial Mass Function (IMF) 
used by us, from lFerrini et alJll992h 
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Fig. 3. a) Dependence of the integrated yield of l4 N produced 
by LIM stars on metallicity for the three yield sets, marked 
with differe nt symbols as in Fig. [2 The crosses are the results 
obtained bv lHenrv et alJ ( 1200 0) for the same VG yields using 
a Salpeter IMF. b) The ratio between the primary to the total 
nitrogen yield with the same symbols than in panel a). The two 
possible techniques to compute the primary component of VG 
yields are represented as short-dashed line (method r) and dot- 
short-dashed line (method AGB). 

While it is ~ 90 % when the AGB technique, represented by a 
dot-short-dashed (blue) line, is used, showing a smaller varia- 
tion with Z. It is interesting that the integrated yield for solar 
abundance in the B U case 3 is around 20%, very similar to the 
value computed by Allo in et al.l 1 11 9791) two decades ago on the 
basis of the observations available at that time. 

All these considerations indicate that the primary nitrogen 
appears at a different time scale in the ISM depending on the 
value of Z. The first primary N will be ejected when stars of 8 
(5 for MA) M Q die, while Z is still low. 

3. The chemical evolution model and its 
calibration 

3.1. Description 

The model used in this wo rk is the Multiphase Ch emical 
Evolution Model des cribed inlFerrini et alJ Jl992i[l994l) . in the 
version presented in Molla & Diaz (2005) and in Paper I. For 
LIM stars, we use the same yields BU than in these two last 
papers, and for comparison purpos es those from MA and VG. 
For massive stars we have chosen Portin ari et alJ y998) and 

3 For Z < 0.0126 we use the same yields set, that is, this one from 
Z = 0.0126. It is always possible to extrapolate the trend obtained for 
the other Z sets, that we show with the straight line joining the points 
corresponding to Z = 0.0126 and Z = 0.0159 



IWooslev & Weaver! i 19951 hereinafter PCB and WW, respec- 
tively). We have run different models computed with different 
combination of yields: BU + WW, VG + WW and MA + PCB, 
which we distinguish as BU, VG and MA, respectively. 

The nitrogen study is usually done by comparing its behav- 
ior relative to iron and oxygen, so it is very important to have 
a careful calibration for these two elements. Oxygen calibra- 
tion was done in Paper I. The SNIa are the main manufacturer 
of iron. The yi elds for type la super nova (SNIa) explosions 
are ta ken from llwamoto et al.l (|1999) and Branch & Nomoto 
(1986). The evolution of this element in the model is quasi- 
independent of the normal stars yields. However, since iron is 
mainly produced by SNIa, even if its yield is very well known, 
its abundance is very dependent on the method to compute the 
rate of these explosions. For this purpose we analyzed the re- 
sults obtained with different possibilities in order to eliminate, 
if possible, uncertainties in the iron abundance evolution. This 
point is relatively important because the Age-metallicity rela- 
tion and the G-dwarf metallicity distribution are usually used 
as calibration methods for chemical evolutions models. In our 
case, furthermore, we compared our results with observed stel- 
lar nitrogen abundances, most of which are given as [N/Fe], and 
so, we checked that the Iron evolution is adequately reproduced 
by our models before this comparison can be made 

We used three methods to compute the SNIa rates 
as given by the following authors: the classical one 
(Matt eucci & Greggiol Il986t iFerrini & Poggiantil Il993h . the 
one given b y iTornambd II1989I) . and other, more recent, 
described in Ruiz-Lapuente et alJ lEoOO). hereinafter named 
MAT, TOR and RL, respectively. The first authors estimate the 
SN rates by using only the Initial Mass Function. The method, 
well described in depth in both cited works, is summarized 
as follows: a proportion of the stellar masses in a given range 
[M„,i„- M max ] will be in binary systems and a fraction of them 
will develop type la supernova. Based on this idea, a mass func- 
tion for the secondary stars is computed from the original one. 
Finally the SNIa rate depends on the number of secondary stars 
that died in each time step, which implies that the time scale for 
the iron appears in the ISM is controlled by the mean lifetimes 
of these secondary stars. 

Actually, this time scale does not depend only on the sec- 
ondary mean lifetimes, since there are other processes that also 
participate in the conversion of a binary system into a SNIa 
explosion. It is necessary to take into account the effects of 
the distances between both stellar components, the orbital ve- 
locities and other parameters to finally obtain the time taken 
for the system to explode since the moment of its formation. 
Tor nambd l ll98 9) performed these calculations for several com- 
binations of possible candidates of binary system or SNIa sce- 
narios (Double Degenerate, Single Degenerate, etc.), provid- 
ing the supernova rate as a function of time normalized for 
a binary system of 1 Mq. All the physical processes and as- 
sumptions are included in their calculations, so we only need 
to include the selected functions in our code and multiply them 
by the number of binary systems, avoiding the need of com- 
puting the secondary and primary initial mass function as de- 
fined in the previously described method. A similar t echnique 
has also been performed more recently bv lRuiz-Lapuente et alJ 
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Fig. 4. Evolution of SNIa/SNII rates for MAT, TOR and RL 
techniques as lab eled in the figure. The ob s erved values for th e 
Galaxy, given by Carjrjellaro et al. ( 1999); Carjpellaro (2004); 
iMannucci et alJ (J2005), are shown, with the error bars, by the 
magenta full dot, the cyan open square and the black triangle, 
respectively. 



(2000). A numerical table has been provided to us by Ruiz- 
Lapuente ( private communication) with the time evolution of 
the supernova rates for a single stellar population, computed 
under updated assumptions about different scenarios and prob- 
abilities of occurrence. We have computed the supernova rates 
using the three methods, thus producing three models MAT, 
TOR and RL. These different techniques affect mostly the iron 
abundances, the other elemental abundances being equal for 
all of them. Therefore, we will compare the three type of SN 
rate calculations by using only BU yields in the analysis of the 
Iron abundance evolution, as well as in the calibration of the 
model (next section). We will compare the three set of LIM 
stars yields when N be studied, using only the RL technique 
and only for the comparison of the relative abundance [N/Fe] 
we will show the nine possible combinations of models. 

The main disparity among the three techniques described 
above resides in the different evolution of the SN rate in 
time. As we see in Fig.|4] MAT is the technique that presents 
highest values of SNIa/SNII at any time, reaching the maxi- 
mum at 2.5 Gyr. RL has a maximum between 2 and 5 Gyr, 
with values approximately 1/2 or 1/3 of those given by MAT. 
Neverthe less it still is within the error bar g i ven bv observa- 
tions Jcappellaro et alJl999tlCappellarol2004llMannucci et alJ 
2005). Note that this value has been reduced for the most recent 
determinations compared with the oldest ones. TOR model is 
the only one with low values. Even if it presents a maximum be- 
fore the first Gyr, this will not be seen in the results because its 
value is very small. From the first Gyr, SNIa/SNII has positive 
slope and it almost reaches the observed value at the present 
time. 



a) 




Time (Gyr) 



Fig. 5. Age-Metallicity relation: a) for BU and RL model di- 
viding WW iron yield by 1, 1.5, 2 and 2.5. b) using MAT, TOR 
and RL SNIa rates with BU yields. Data are from authors of 
Table^as labeled. 

3.2. Calibration of the model: Iron evolution in the 
Solar Vicinity 

The results for iron abundance obtained with these three meth- 
ods are shown in Fig. [5] the Age-Metallicity Relation (AMR) 
and in Fig.|6j the G-Dwarf distribution, for the Solar Vicinity. 
For this comparison we have shown only BU yields, keeping 
in mind that the set of yields will have only small effects on 
this relation. Nevertheless, this model uses WW yields for mas- 
sive stars and these authors claimed that th is cou ld produce too 
much iron, and advised in Tim mes et alJ i 19951) to divide the 
iron ejections at least by two. In order to calculate how much 
that WW iron excess will be, we have calculated four differ- 
ent models for BU yields and RL technique, where massive 
stars iron production is divided by 1, 1.5, 2 and 3. Results are 
represented in panel a) where it is clearly shown that a fac- 
tor of 2 is a good compromise that we will use in panel b). In 
this last panel we present the Age-Metallicity relation for the 
three SNIa cases, where all of them are in reasonable agree- 
ment with data, given their large dispersion. Although there are 
very little differences between models, it can be seen that the 
iron appears later and takes a little more time to reach high val- 
ues when MAT and RL techniques are used, than for the TOR 
SNIa method, but all of them reach the solar abundance. 

Regarding the G-dwarf distribution, represented in Fig.|6] 
differences appear mainly between TOR and the others because 
it provides a narrower distribution than the others. The three 
models are able to reproduce the low metallicity tail without 
showing any G-dwarf problem. 

In Fig. we show the relation between iron and oxygen. 
As before, in panel a) the BU + RL model is presented varying 
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X 


— 


— 


Reddv et al. (2003)(RED) 


X 


X 


X 


— 


— 


X 


Rocha-Pinto et al. (2000b.a)(RO) 


X 


— 


— 


— 


— 


X 


Rolleston et al. (2000), 


X 


X 


X 


X 


X 


— 


Smartt & Rolleston (1997). 














Smartt et al. (2001) 














Shi et al. (2002) 


X 


X 


X 








Smith etal. (2001) 


X 






X 






Smte et al. (2005)(SPI) 


X 


X 


X 


X 






Tomkin & Lambert (1984). 


X 


X 


X 


X 






Tomkin et al. (1986. 1995). 














Twarog (1980)(TW)~ 


X 










X 


Westin et al. (200O) 


X 


X 


X 


X 







Table 1. References for CNO stellar abundances used for the comparison with model results. 



the massive stars iron ejection. In this case the differences are 
clearer than in the AMR case. We chose the model Fe/2 that 
we will use for the rest of the paper. In panel b) we plot the 
model results using BU yields with the three SNIa techniques. 
As in the previous case, the LIM stars yields do not change the 
results because oxygen is ejected by massive stars and iron is 
mainly produced by SNIa events. The three models have very 
similar behavior. 



4. Results analysis: The nitrogen abundances 

We devote this section to analyzing the results obtained with 
our models for nitrogen abundances and comparing them with 
the available observational data. We divide these results in four 
parts: a) the evolution of Nitrogen over time for the Solar 
Vicinity (assumed located at a galactocentric distance of 8 kpc), 
b) the radial distributions of elements in the disc, c) the relation 



of log(N/0) with the oxygen evolution and d) the relation of 
[N/Fe] with iron. 

4.1. Time evolution of nitrogen 

Fig.[8]shows the evolution of nitrogen with t he solar abundance 
values -large fil led symbols- tak en from Grevesse & Sauval 
( 1998) - circle-, iHolwegeilfcOOll) -square- and Asnlun det alJ 
( 2005), -cross-, by assuming an age of 4.5 Gyr for the sun. For 
the interstellar medium abundances at 13. 2 Gyr, large empty 
symb ols, we use th e abundances giv en by Mev er et all £1997, 
1998), -circle-, and Peimbert (1999) -square-. The small open 
circles are abundances for objects with given stellar ages, at a 
radial galactocentric distance between 7.5 and 9.5 kpc. 

Models BU, VG and MA are represented by the (red) solid, 
the (blue) short-dashed and the (green) long-dashed lines, re- 
spectively. Both models VG, following the two possible tech- 
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Fig. 6. T he G-Dwa rf distribution in the Solar Vicinity. Data ar e 
from Chang et al. (2000) (dashed line) and J0rgensen (2000), 
(long-dashed line) 



niques to calculate the proportion of primary nitrogen, tech- 
niques r and AGB described in Section |2] yield results indistin- 
guishable for times larger than 1 Gyr, so we represent only the 
results for the first one. In panel a) MA and VG models give a 
greater value than BU since at the lowest metallicity their nitro- 
gen yield are higher than the corresponding one from BU (see 
FigGJ. Then, once an abundance higher than ~0.004 is reached, 
it continues increasing smoothly until the present time, repro- 
ducing both, solar and ISM, abundances. The shape shown by 
the three models are similar and all of them reproduce both the 
solar value and the ISM value. 

The same kind of information can also be extracted from 
the relative abundances represented in panel b). In panel b), 
we show the time evolution of log(N/0). Since there is a good 
agre ement in fitting the abundance of oxygen for all models, 
(see iGavilan et aljl2005t) . the differences in this plot must be 
only due to the nitrogen production. The disagreement between 
the different models is important for times shorter than 1 .5 Gyr, 
when intermediate mass and massive stars are the main con- 
tributors and the distinct primary/secondary ratios effects are 
evident there. Model MA has a strong increase in the first Gyr 
due the primary component and then it flattens up. Model BU 
has a ratio of primary nitrogen larger than MA for all Z except 
for the lowest one, what produces a smoother evolution, and, 
finally, N remains below MA. The resulting final N/O ratios 
are similar in both models and in agreement with observations. 
The behavior mostly primary of model VG, when all AGB ni- 
trogen is considered as such, implies a very strong increase of 
the abundance at the earliest times. After that, both methods 
give a smooth slope, reaching an absolute value around -0.5 
dex higher than the observed one. 



Fig. 7. The relation [O/Fe] vs [Fe/H] : a) The model results for 
BU + RL varying WW iron ejection; b) model BU using MAT, 
TOR and RL SNIa techniques and BU yields. Open (black) 
dots are stellar data from authors listed in Table [J from w hich 
we re presented the mo st recent from Israe lian et al.l (12004) and 
lAkerman et alJ J2004I) as (red) triangles, and full (blue) dots. 



The good behavior of BU yields is also evident in panel c) 
where log(N/C) is shown. MA model presents a maximum at 
the first Gyr, the decrease is due to the higher amount of carbon 
ejected in that model (see Paper I); so the absolute value at 
the present time is only marginally reached. The shape for VG 
model is similar than the BU one but with a nitrogen excess. 
All models seem to fit the solar and ISM data but the model 
BU is the best one in reproducing the stellar data, and, more 
importantly, it is the best one at fitting all of the data at the 
same time. 

4.2. Nitrogen abundance in the Galactic disc 

Now, we will explore the radial distributions of nitrogen over 
the galactic disc, shown in Fig. [9] Data correspond to Hn re- 
gions from references as labeled in the figure and to stars from 
references in Tabled 

The radial distribution is more or less reproduced within 
the errors by all models. Actually, the shape of the radial dis- 
tribution is well fitted in all cases, independently of the abso- 
lute values, since this results as an effect of the ratio infall/SFR 
along the galactocentric radius, produced by the scenario of our 
MWG model, and therefore is rather independent of the yields 
used. However, the slope of the radial distribution at the two 
ends of the disc, in the c enter and in the oute r region s, is a mat- 
ter of discussion. Thus, IVilchez & Estebanl 1 1996) claim that 
the gradients are not as steep in these regions as in the rest of 
the galactic disc. The same occurs in the inner disc where the 
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Fig. 8. Time evolution of elemental abundances in Solar 
Vicinity for a) nitrogen, as l2+log{N/H);b) nitrogen over oxy- 
gen as log(N/0) and c) nitrogen over carbon as log(N/C). Solar 
abundances are the filled symbo ls from iGrevesse & Sauvall 
1 19981) - ckcle-. lHolwe"gerlll200ll) -square- and Asplund et al. 
(2005), -cross-, by assuming an age of 4.5 Gyr for the 
sun. Empty symbols at 13.2 Gyr are the inters tellar medium 
abundanc es give n by iMever et alJ i 19971 Il998l) . -circle- and 
IPeimberj dl999l) -square-. Small open dots are stellar abun- 
dances obtained by authors from Table [2 being those located 
around the Solar Vicinity (7.5 < R < 9.5 kpc). Lines meaning 
is given in panel a). 



Fig. 9. Radial distributions of nitrogen abundances, a) as 12 
+ log (N/H), b) as log(N/0) and c) as log(N/C) for the 
same models than the previous figure. Stellar data -open 
dots- are taken from author s given in Ta ble [H including the 
most recent ones from iDaflon & Cunhal (2004) marked sep- 
arately. G alactic Hn regio ns ab undances are also included, 
taken fromlPeimberll i ll 9791) -PE-.IShaver et alldl983h -SHA- 
iFich & Silkevl(ll994 -FS- . IVilchez & Estebanl ill 9961) -TOS- 
lAfflerbach etalJ (Il997l) -AFF-. and lEsteban etal] i 1 9991 
2005l): ICarigi et alJ ( 120051) -CAR- as marked in the figure. The 
solar and interstellar abundances have the same symbol than in 
Fig. El 



most recent data from Smart tetafl J200ll) show that the distri- 
bution flattens. Our models have been tuned to fit these two sets 
of data, and due to that, the resulting overall gradient is smaller 
than that obtained by other authors. As can be seen in Fig. 
MA and VG models produce a gradient flatter than observed. 

In Fig. |9j>) and c), the radial distributions log(N/0) and 
log(N/C) are plotted, as they are considered important for the 
study of different yields. In panel b) the radial distribution of 
log(N/0) showed by data presents a clear slope, although there 
is some data that shows a flatter distribution in the outer re- 
gions. A steep radial distribution for N/O is expected because 
oxygen is produced by massive stars. If nitrogen was ejected 
by massive stars, its secondary character would cause it to en- 
ter the ISM after the oxygen. Instead, if it would be ejected by 
intermediate stars, the time needed for their evolution would 
be larger. Thus, in both cases, the nitrogen appears in the ISM 
after the oxygen does. 

Once again, the large dispersion of the data prevents a 
clear selection of the best model, however the BU model 
seems most adequate to r epro duce the Hn regions data from 
IVilchez & Estebanl i 19961) and lFich & Silkevl fl99ll) . The MA 
model shows a radial gradient flatter than indicated by observa- 
tions, with higher absolute values compared to the mean values 



of data, mostly in the outer disc, and VG models, as before, 
show higher values. The same arguments are also valid when 
panel c) is analyzed. In this case a slight negative gradient is 
shown for log(N/C). The small amount of data and the great 
dispersion prevent the selection of any model as better than the 
others, although it seems clear that MA remains below most of 
them as corresponds to the large production of C, and that VG 
lies in the upper side of the data. It is apparent that BU model 
shows a better behavior compared to the data. Once again we 
stress the importance of using adequate yields to reproduce the 
whole set of data at the same time. Yield BU seems to be in the 
adequate range of production of N, C and O, since the model 
appears in the zone occupied by of data in the three panels. 

It is necessary to remember that open dots represent stellar 
abundances. We have tried to select only those corresponding 
to young stars, but we do not know the age of the complete set 
of stars with available data. In this case we have preferred to 
use the available abundances; thus, it is possible that some data 
does not correspond to young enough stars. 
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4.3. Nitrogen vs Iron 

In this case, as the iron evolution may have also an influence 
over the model results, we have presented the relation between 
nitrogen and iron, Fig.^3 with a different panel, a), b) and c), 
for each set of yield, MA, BU and VG, respectively. The three 
possible methods to compute the SNIa, RL, TOR, and MAT are 
shown with solid (red), short-dashed (blue) and long-dashed 
(green) lines, respectively, in each panel. The first thing we ob- 
serve is that the effects of the different SNIa techniques are 
almost indistinguishable. Therefore, the main features of each 
model at those metallicities are due to yields. In other words, 
we may analyze the behavior of the nitrogen corresponding to 
each yield set disregarding the accuracy in the SNIa calcula- 
tions. 

All results agree in the sense that the first nitrogen to be 
ejected is secondary, as due to the massive stars, so the ini- 
tial slope is positive and large, -although this behavior is not 
shown in the figure because it occurs when oxygen abundances 
is lower that 5 dex- but they differ in when the slope begins 
to change. When the N ejected by LIM stars appear, there is 
a strong increase due to the change from a secondary to a pri- 
mary behavior. In MA yields LIM stars eject less primary ni- 
trogen, and later, since it is ejected as secondary for stars up 
to 5 M Q . Then the main contributors to primary N are the stars 
with masses between 2M© and 3M©. For this reason the slope 
does not change until it reaches [Fe/H] = -1 .5, the moment in 
which these stars begin to die. When BU yields are used, the 
trend changes earlier in the evolution due to the contribution 
of the primary nitrogen ejected by stars in the range 4-8 M . 
As their life is so brief, the ejection occurs at [Fe/H] = -4. 
From then, the slope is close to zero: the signature of primary 
nitrogen. The case of VG shows a behavior similar to BU. 

We would like to remark that the data dispersion is so great 
that all the models lie in the data area, regardless of their big 
discrepancies, although the region of the metal-rich objects 
{[Fe/H] > -1.5) is particularly well fitted in panel b) by Model 
BU. It is necessary to use the very low metallicity data to clar- 
ify which mo del works better. In f act, the most rece nt obser- 
vations from llsraelian et all d2004t) : ISoite et all d2005l) show a 
slope flatter than before, (even with a negative slope) which 
is a behavior more consistent with our model BU than with 
the model obtained with MA yields. This last model might be 
considered acceptabl e when the available low metallicity data 
were only those from Carb orTet alJ dl987h . but when using the 
new determinations of N abundances for this kind of objects, 
the conclusion is that BU better reproduces the generic trend of 
data. It is also necessary to take into account that most of the 
metal-poor objects do not belong to the disc but to the halo. In 
this way, we represent the halo model results for the zone that 
infall over the disc at galactocentric distance equal to 8 Kpc, at 
the right panels b), d) and f) of Fig^H W e see that the tren d 
shown by t he rece nt observations from Israelia n et al 1 d2004l) : 
Spit e et all d2005l) is more compatible with BU and VG than 
MA. 



, , , , 1 , , , , 1 , , , , 1 , , , , r 

- BU 

MA 

VG-r 

- VG-AGB 




6 7 8 9 

12 + log(0/H) 



Fig. 11. The relative abundances of log (N/O) vs the Oxygen 
abundance as 12 + log(0/H) for the Solar Neighborhood (lines 
as shown in the plot). Stellar abundances are the open dots 
from Table HI from whi ch we mark sep arately those from 
Israelian et al. (2004) and Spi te et al] {2005) as (red and ma- 
genta) full triangles. Galactic Hn regions data in th e Sola r 
Region (7.5 < R ( kpc) < 9.5 are from IPeimberd d 19791) 
IShaver et alJ dl983UFich & Silkevl |l99l.lVilchez & EstebarJ 
dl996l) . iTsamis et alJ d2003l) and Icarigi et alJ d2005l) . The 
meaning of symbols is given in the figure. The solar position 
is marked with arrows. 

4.4. Nitrogen vs oxygen 

Finally, we show in Fig. the classical and well known 
graphic of the relative abundance of N vs O as log(N/0) 
vs 12 + log(0/H). We show the final results for the Solar 
Neighborhood of the computed models. Model BU reproduces 
well the expected behavior of N when the whole data is taken 
into account. Not only the level of N is so adequate but also the 
shape is smoother than the one shown by the other two models. 
The observed trend at low metallicity can be well reproduced 
with BU yields because they have the appropriate primary to 
secondary ratio, and the adequate integrated nitrogen yields. 
MA yields have also a primary nitrogen component, but the in- 
tegrated nitrogen yield has a metallicity dependence in the op- 
posite way as BU for the lowest Z, so the trend shown by data 
can not be well reproduced. Both VG models have the right 
shape and are almost the same for 12 + log (O/H) > 8. The 
problem is that the integrated yields are high. It would be nec- 
essary to change the input parameters as the infall rate or the 
efficiencies to form stars in order to fit the solar abundances. In 
that case, probably, other data will not be reproduced. Only the 
BU model shows simultaneously the good shape and adequate 
absolute abundances. For comparison purposes, we have also 
shown the resulting model using VG yields but assuming that 
the Nitrogen is completely secondary. 

If the flat behavior were caused by massive stars, the data 
dispersion would b e very small. A problem arises when the 
metal-poor objects ( Israelia rTet alJl2004l ISpite et aTll2005l) are 
included in the figure, as can be seen in Fig. \H\ Some values 
follow the described trend over the flat line, but there exist some 
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Fig. 10. The relative abundances of [N/Fe] vs the iron abundance [Fe/H] for the solar region. The panels on the left show the disc 
evolution and the panels on the right show the halo evolution. In each one of them, the evolution obtained with three different 
methods to compute the SN-Ia rates are show n as labeled in panel a ) . Data are taken from the different referenc es listed in Tabled 
marked separately the most recent ones from llsraelian et al T j2004 : lEcuvillon et alJ ( Eooi : ISpite et alJ ( 120051) in panel d). 



lower abundances which are around log(N/0) 2. This be- 
havior is not compatible with a primary component proceeding 
only from massive stars. 

We show in Fig^l me evolution given by BU model for 
four different radial regions of the Galaxy: two inner (~ 2 and 4 
kpc) more evolved regions -(magenta) dotted and (blue) short- 
dashed lines-, the solar vicinity (~ 8 kpc) as in the previous 
figure, the (red) solid line, and an outer one (~ 18k pc), - 
the (green) long-dashed line-, where the evolution takes place 



slowly. In panel a) we show the results for the halo zones and in 
panel b) for the disk regions. In both panels we have included 
the stellar data for the MWG. 

The numbers on the graph indicate the evolutionary time, 
in million years, that corresponds to that point of the line. This 
is necessary because the 12 + log(0/H) value is not the same 
for each radius at the same value of time. 

The halo regions have similar evolutions independent of 
their distance from the center of the Galaxy. All of them 
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R=2 kpc 




12 + log(0/H) 



Fig. 12. The relative abundances of log (N/O) vs the Oxygen 
abundance as 12 + logiO/H) for the halo -upper panel- and 
for the disc -lower panel- for four different galactocentric re- 
gions: two inner regions located at 2 and 4 kpc, dotted and short 
dashed lines, a solar region, solid line, and an outer region, long 
dashed line, respectively. Data refer to Galactic stars are from 
the same authors as cited in Tabled The meaning of symbols 
is given in the figure. The solar position is marked with arrows. 

reproduce well t he recent data from Israelia rTet alJ (|2004); 
Scit e et alJ J2005I) obtained for 5.5 < 12 + log(0/H) < 8, and 
the disk regions fit the stellar data obtained for 12 + log{0/H) > 
8 of the disk. Their evolutionary tracks, however, are very dif- 
ferent, as corresponds to their distinct input parameters (infall 
rates, initial gas masses, efficiencies to form stars...) which are 
translated into very different star formation histories. Thus, the 
dispersion of the MWG data can be well explained on the ba- 
sis of a primary production of nitrogen from LIM stars, higher 
for the lowest metallicities, and with different star formation 
efficiencies in the different regions. 

We represent the same results for the disk regions in 
FigE] compared with data referring to Galactic HII regions, 
taken from the same authors than those of Fig. but with- 
out li miting the possible g a lactoc e ntric distance. Other g alaxy 
data JGarnett et alJ 1 19951 Il999t Ivan Zee et al] |l998a) and 
llzotov&Thuanl (ll999) are also shown. The large open triangle 
is the recent estimate obtained from llzotov et alJ d2005l) for the 
lowest metallicity known galaxy. We have also added the DLA 
objects data fromlpettini et all J2002l) : IProchaska et all 0002): 
ICenturion et alJ J2003) as solid points. We want to remark that 
the disk regions evolve in good agreement with all of them, 
showing the inner regions a steeper evolution while the outer 
one shows a very fiat evolution with a high and constant value 
log(N/0) 1.2 dex, similar to the behavior of dwarf galax- 
ies. These results suggest that the observed dispersion in this 
kind of plot, when other galaxies data (such as dwarf or DLA 
galaxies) are included, might be reproduced if different star for- 
mation histories have occurred in different galaxies. This argu- 
ment Jias_ajready_been in voked by other a uthors, in particular 
bv lHenrv et ailfcOOOl) andlPrantzosI J2003b . It was even demon- 
strated bv lPilvugin et alJ d2003l) . who analyzed data for differ- 




[Fe/H] 



Fig. 13. The relative abundances of log (N/O) vs the Oxygen 
abundance as 12 + log(0/H) for the same disk regions than 
previous figure compared with data referring to Galactic and 
extragalactic HII regions and DLA objects data. The meaning 
of symbols is given and the solar position is marked with ar- 
rows. 

ent radial regions in spiral galaxies and showed the changes of 
the evolutionary track in the plane N/O vs O/H for each one of 
them. It is evident that this kind of behavior may be represented 
by our models and that the new yields may reproduce better the 
whole set of data. In fact, these yields have already been used 
in a grid of chemical evolution models for a large number of 
theoretical galaxies [Molla & Dfazl J2005I) . A discussion about 
the resulting N/O abunda nces and its possib le dispersion for 
different objects is done in Molla et all J2005t) . 

This figure and the behavior of (C/O) vs O/H, shown in 
Paper I, are the main clues to consider the present yields the 
most adequate to represent the evolution of galaxies. The pro- 
duction of carbon by LIM stars is sufficient to obtain an in- 
crease in C/O without the need to invoke mass loss by massive 
star winds, and the N/O behavior may be well reproduced with 
different star formation efficiencies due to the adequate level of 
the primary component produced by LIM stars and to the right 
dependence of this component with Z. 

5. Conclusions 

Our conclusion can be summarized as follows: 

1 . The primary component of nitrogen, necessary to explain 
the trend of N/O with O/H, may be mostly produced by 
LIM stars, and adequately fits all the data, including the ob- 
served dispersion. In this way the integrated yield produced 
by LIM stars must be directly proportional to Z, while the 
ratio NplN tot must increase for Z decreasing. A primary 
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component, larger for lowest metallicities, has an important 
effect on explaining low abundance range data. 

2. The dependence of the N yield on stellar mass would have 
a maximum around 5-6 M Q , while the primary compo- 
nent shows other around 3.5-5 M . This constraints the 
time where these contributions have important effects on 
the evolution. 

3. The high dispersion on N/O data for low and high metal - 
licity galactic regions may be explained with these yields 
as we have demonstrated with the radial regions of the 
disc models which have different star formation efficien- 
cies. O ur fi ndings go in the same address than lHenrv et alJ 
(2000) and Prantzos (2003) using VG yields, but we claim 
that it is easier to reproduce the whole data set when BU 
yields are used. Our model for MWG with small efficien- 
cies to form stars is consistent with data. 

4. These results seem suggest that models with differences 
in the star formation histories for diff erent types of galax- 
ies, such as those calculated in Molla & Diaz (2005), might 
produce final abundances with high dispersion, in agree- 
ment with the observed one when dwarfs galaxies or DLA 
galaxies are included in a plot N/O-O/H, such as we will 
sh ow inlMolla et alJd2005l) 

5. AslChiaPDin i et al.l (120 03b). we also support that the halo 
and the disc have different evolutions. The set of stellar data 
[C/Fe], [N/Fe] and [N/C] may be divided into two trends. 
The first one is well reproduced by our disc models, while 
the second one is well fitted by our halo results. 

6. Due to the primary N component of BU yields, and since 
the intermediate stars have short lifetimes, it is possible 
to produce high [N/Fe] abundances even at low metallic- 
ities which are in per fect agreement with the recent halo 
stars data o btained by Israeli an et al 1 d2004l) : ErvTeT et al. 
j2004l) : LSpite et alJd2005l) . 

7. We claim that the fit of the whole set of data with only one 
model is not an easy task. We may reproduce the observed 
trend with BU yields combined with yields from WW. In 
summary, our model BU reproduce reasonably well the 
whole CNO data set. 
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